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The results ol a search for squarks and gluinos using data from pp collisions recorded at a center- 
of-mass energy of 1.96 TeV by the D0 detector at the Fermilab Tevatron CoUider are reported. 
The topologies analyzed consist of acoplanar-jet and multijet events with large missing transverse 



energy. No evidence for the production of squarks or gluinos was found in a data sample of 310 pb~^. 
Lower limits of 325 and 241 GeV were derived at the 95% C.L. on the squark and gluino masses, 
respectively, within the framework of minimal supergravity with tan/? — 3, Aq — 0, and /i < 0. 

PACS numbers: 14.80.Ly, 12.60.Jv, 13.85.Rm 



Supersymmetric models predict the existence of spin-0 
quarks, or squarks (g), and spin-1/2 gluons, or gluinos 
(^), as partners of the ordinary quarks and gluons. Su- 
persymmetric particles carry a value of —1 for i?-parity, 
a multiplicative quantum number, while R = 1 for stan- 
dard model (SM) particles. If i?-parity is conserved, 
as assumed in the following, supersymmetric particles 
are produced in pairs. Their decay leads to SM parti- 
cles and to the lightest supersymmetric particle (LSP), 
which is stable. In supersymmetric models inspired by 
supergravity ^^, the commonly accepted LSP candidate 
is the lightest neutralino {xi, a mixture of the super- 
partners of the neutral gauge and Higgs bosons), which 
is weakly interacting, thus escaping detection and pro- 
viding the classic missing transverse energy (^t) signa- 
ture at colliders. The most copiously produced super- 
symmetric particles in pp collisions should be, if suffi- 
ciently light, colored particles, i.e. squarks and gluinos. 
If squarks are lighter than gluinos, they will tend to de- 
cay according to g — > gxi? ^^d their pair production will 
yield an acoplanar-jet topology wither- If gluinos are 
lighter than squarks, their pair production and decay via 
g — » qqxi will lead to topologies containing a large num- 
ber of jets and^T- 

In this Letter, a search for squarks and gluinos in 
topologies with jets and large I^t is reported, using 
310 pb~^ of data collected at a center-of-mass energy of 
1.96 TeV with the D0 detector during Run II of the Fer- 
milab Tevatron pp Collider. The search was conducted 
within the framework of the minimal supergravity model 
(mSUGRA) [l|. Previous direct mass hmits are 195 GeV 
for gluinos if squarks are very heavy, and 300 GeV for 
squarks and gluinos of equal masses |3, yf ■ 

A detailed description of the D0 detector can be found 
in Ref . [J| . The central tracking system consists of a sil- 
icon microstrip tracker and a central fiber tracker, both 
located within a 2 T superconducting solenoidal magnet. 
A liquid-argon and uranium calorimeter covers pseudora- 
piditics up to \ri\ « 4.2, where 77 = — In [tan (0/2)] and 9 
is the polar angle with respect to the proton beam direc- 
tion. The calorimeter consists of three sections, housed 
in separate cryostats: the central one covers \ri\ < 1.1, 
and the two end sections extend the coverage to larger 
?7|. The calorimeter is segmented in depth, with four 
electromagnetic layers followed by up to five hadronic 
layers. It is also segmented in projective towers of size 
0.1 X 0.1 in ?7 — space, where (f) is the azimuthal an- 
gle in radians. Calorimeter cells are defined as intersec- 
tions of towers and layers. Additional sampling is pro- 
vided by scintillating tiles in the regions at the bound- 



ary between cryostats. An outer muon system, cover- 
ing \r]\ < 2, consists of a layer of tracking detectors and 
scintillation trigger counters in front of 1.8 T toroids, 
followed by two similar layers after the toroids. Jets 
are reconstructed from the energy deposited in calorime- 
ter towers using the Run II cone algorithm [5| with ra- 
dius 7?. = ■y/(A0)2 -I- (A?7)2 =0.5. The jet energy scale 
(JES) is derived from the transverse momentum balance 
in photon-plus-jet events. The ]^t is calculated from all 
calorimeter cells, and corrected for the jet energy scale 
and for reconstructed muons. 

The D0 trigger system consists of three levels, LI, 
L2, and L3. The events used in this analysis were 
recorded using a jet trigger requiring missing transverse 
energy calculated using the sum of the jet momenta 
{]^T = I J2\ctsPT\)- At LI, events were required to have at 
least three calorimeter towers of size A(/) X A77 = 0.2x0.2 
with transverse energy Et greater than 5 GeV. Events 
with a large imbalance in transverse momentum were 
then selected by requiring ^t to be greater than 20 GeV 
and 30 GeV at L2 and L3 respectively. In a small fraction 
of the data sample recorded at a higher instantaneous lu- 
minosity, the acoplanarity, defined as the azimuthal angle 
between the two leading jets, was required to be less than 
168.75° and 170° at L2 and L3 respectively. 

The signal consists of jets and ^t- This topology 
also arises from SM processes with real ]^t, such as 
pp — > Z-f jets with Z — > vD, and from multijet pro- 
duction when one or more jets are mismeasured (QCD 
background). Simulated events from SM and mSUGRA 
processes were produced using Monte Carlo (MC) gen- 
erators, subjected to a full GEANT-based jg simulation 
of the detector geometry and response, and processed 
through the same reconstruction chain as the data. The 
CTEQ5L [3] parton density functions (PDF) were used, 
and a Poisson-average of 0.8 minimum bias events was 
overlaid on each simulated event. The QCD background 
was not simulated, but estimated directly from data. To 
simulate WVZ-l-jets and ti production, the ALPGEN 1.3 
generator Q was used, interfaced with PYTHIA 6.202 [3 
for the simulation of initial and final state radiation and 
of jet hadronization. The next-to-leading order (NLO) 
cross sections were computed with MCFM 3.4.4 [lO|, or 
taken from Ref. [ll| for ti production. 

Squark and gluino production and decay were simu- 
lated with PYTHIA. The masses and couplings of the su- 
persymmetric particles were calculated with is A jet 7.58 
J12l | from the set of five mSUGRA parameters: ttiq and 
TOi/2, which are universal scalar and gaugino masses, and 
Aq, a universal trilinear coupling, all defined at the scale 



of grand unification; tan /3, the ratio of the vacuum ex- 
pectation values of the two Higgs fields; and the sign of 
the Higgs-mixing mass parameter /i. To retain consis- 
tency with earlier analyses [^ |3|, the following param- 
eters were fixed: ^o = 0, tan/3 = 3, and ^ < 0. For 
the same reason, the production of scalar top quarks, 
or stops, was ignored. In the following, "squark mass" 
stands for the average mass of all squarks other than 
stops. All squark and gluino decay modes were taken 
into account in the simulation, including cascade decays 
such as g ^ qqX2 with X2 ~^ ^~^('~Xi- The NLO cross 
sections of the various signal processes were calculated 
with PROSPINO 2 [lj|. 

Three benchmark scenarios have been considered. At 
low mo, the gluino is heavier than the squarks, and the 
process with the dominant cross section is qq production. 
A "dijet" analysis was optimized to search for events con- 
taining a pair of acoplanar jets. At high mo, the squarks 
are much heavier than the gluino, and the process with 
the highest cross section is therefore gg production. A 
"gluino" analysis was optimized to search for multijet 
events (> 4 jets). In the intermediate mo region, all 
squark-gluino production processes contribute to the to- 
tal cross section, in particular the qg process becomes 
relevant. A "3-jets" analysis was optimized to search for 
events with at least three jets. The benchmark for this 
analysis is the case where rriq — rrig. 

A common event preselection was used for the three 
analyses to select events with at least two jets and sub- 
stantial^T (> 40 GeV). The acoplanarity was required to 
be below 165°. The longitudinal position of the primary 
vertex with respect to the detector center was restricted, 
\z\ < 60cm, to ensure an efficient primary vertex recon- 
struction. The two leading jets, i.e. those with the largest 
transverse energies, were required to be in the central re- 
gion of the calorimeter, |?ydct| < 0.8, where rydct is the 
jet pseudorapidity calculated under the assumption that 
the jet originates from the detector center. These jets 
must have their fraction of energy in the electromagnetic 
layers of the calorimeter smaller than 0.95. Minimum 
transverse energies of 60 and 40 GeV were required for 
the first and second leading jets, respectively. 

The tracking capabilities of the Run II D0 detector 
were used to significantly reduce the QCD background. 
A comparison of the jet energy with the energy carried 
by its associated charged particles was performed. In 
particular, the ratio CPF of the transverse momentum 
carried by tracks associated with the jet to the jet Et 
is expected to be close to zero if an incorrect primary 
vertex was selected. The two leading jets were required 
to have CPF larger than 0.05. 

Different selection criteria were next applied in the 
three analyses, as summarized in Table Q] In the "di- 
jet" analysis, the cut on the second jet Et was raised to 
50 GeV. In the "3-jets" and "gluino" analyses, a third and 
fourth jet were required, respectively. They must fulfill 



TABLE I: Selection criteria for the three analyses (all energies 
in GeV); see the text for further details. 



Preselection Cut 




All Analyses 




Acoplanarity 
[Vertex z pos.j 




>40 
< 165° 
< 60 cm 




Selection Cut 


"dijet" 


"3-jets" 


"gluino" 


1st jet Et" 
2nd jet Et'" 
3rd ]et Et" 
4th jet Et" 


>60 
>50 


>60 
>40 
>30 


>60 
>40 
>30 
> 20 


Electron veto 
Muon veto 


yes 
yes 


yes 
yes 


yes 
yes 


A0(^T,jeti) 

A(?i(^T,jet2) 

A(?!)mm(^T,anyjet) 


>90° 
> 50° 
>40° 


>90° 
>50° 


>90° 
> 50° 


Ht 

IpT 


> 275 

> 175 


> 350 

> 100 


> 225 
> 75 



"Jets subject to an Et cut are also required to be central 
(I'JdctI < 0-8)i with an electromagnetic fraction below 0.95, and to 
have CPF>0.05. 



the same quality criteria as the two leading jets, except 
for the Et cuts which were set at 30 and 20 GeV. In all 
three analyses, a veto on isolated electrons or muons with 
PT > 10 GeV rejects a large fraction of events originating 
from the W/Z+]eis processes. The azimuthal angles be- 
tween the^T and the first jet, A(/)(^r, jeti), and the sec- 
ond jet, A0(^7-, jet2), were used to remove events where 
the energy of one jet was mismeasured, generating IpT 
aligned with that jet. The cuts are A(^(^7-,jeti) > 90° 
and A0(^T,jet2) > 50°. 

In the "dijet" analysis, QCD events were further sup- 
pressed by requiring that the minimum azimuthal an- 
gle A0mi,i(^T, anyjet) between the^^T and any jet with 
Et > 15 GeV be greater than 40°. Because of the higher 
jet multiplicity, this criterion was not used in the "3-jets" 
and "gluino" analyses. 

The "dijet" A0niin(^T, anyjet) cut along with the two 
final cuts on Ht — X^icts-^T and on^r were optimized 
by minimizing the expected upper limit on the cross sec- 
tion in the absence of signal. To this end, as well as for 
the derivation of the final results, the modified frequentist 
CLg method |lj| was used. For each set of cuts tested, 
the QCD background contribution was estimated from 
an exponential fit to the ^t distribution below 60 GeV, 
after subtraction of the SM background processes, ex- 
trapolated above the chosen ^t cut value. The optimal 
cuts thus determined are given in Table Q] for the three 
analyses. Figure Q] shows: the A(/)inin(^T,anyjet) distri- 
bution after applying the "dijet" analysis criteria with a 
IpT cut reduced to 80 GeV and without requiring the con- 
ditions on A(/)inin(^T,anyjet) itself and on A(/)(^T,jet2); 
the Ht distribution after applying all the "3-jets" analy- 
sis criteria except the one on Ht] and the IpT distribution 
after applying all the "gluino" analysis criteria except the 
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FIG. 1: A<^niin(^T, anyjet) distribution after applying the "dijet" analysis criteria with a _^t cut reduced to 80 GeV and 
without requiring the conditions on A(/)inin(^T,anyjet) itself and on A(?i(^T, jet2) (left), Ht distribution after applying all the 
"3-jets" analysis criteria except the one on Ht (middle), and.^T distribution after applying all the "gluino" analysis criteria 
except the final one on ]^t (right), for data (points with error bars), for non-QCD SM background (full histogram), and for 
signal MC (dashed histogram on top of SM). For each analysis, the signal drawn is the one for the appropriate benchmark 
scenario (Table ITB. In thc^r distribution, the fitted QCD background is also drawn. 

TABLE II: For each analysis, information on the signal for which it was optimized: mo, mi/2, rrig, rriq and nominal NLO cross 
section, signal efficiency, the number of events observed, the number of events expected from SM and QCD backgrounds and 
the 95% C.L. signal cross section upper limit. The first uncertainty is statistical and the second is systematic. 



Analysis 


(mo,mi/2) 
(GeV) 


(mg,mq) 
(GeV) 


(Pb) 


(%) 


A^'obs. 


Abackgrd. 


0-95 

(pb) 


"dijet" 

"3-jets" 

"gluino" 


(25,145) 
(191,126) 
(500,80) 


(366,318) 
(330,330) 
(240,507) 


0.63 
0.64 
2.41 


&.2±Q.4+ll 
4.7 ± 0.3t°;« 
2.3±0.2t;5;3 


6 
4 
10 


■ 0+4.4 +1.1 
^■°-2.0 -0.8 
q Q+1.3 +0.7 
^■^-1.0 -0.8 

in ■?+i-5 +1.9 

iU.d_l 4 _2.5 


0.44 
0.45 
1.72 



one on^T- 

The numbers of events selected by each analysis are re- 
ported in Table ITU as well as the numbers of background 
events expected. Six events were selected by the "dijet" 
analysis, four by the "3-jets" analysis, and ten by the 
"gluino" analysis. The total expected background con- 
tributions are 4.8, 3.9 and 10.3 events, respectively. The 
main background contributions are from Z — > i/f' -I- jets, 
W ^f Iv -I- jets, and ti^ hh qq'lv. The QCD background 
was evaluated from a fit to the Ipr distribution as de- 
scribed above. It was found to be negligible in the "dijet" 
and "3-jets" analyses, and was therefore conservatively 
ignored. A QCD contribution of 0.7lo;4 event was esti- 
mated in the "gluino" analysis. The uncertainties were 
obtained by taking into account the accuracy of the fit 
parameter determination and by varying the range of the 
fit. The signal efficiencies are given in Table HTl for the 
three benchmark scenarios, with the corresponding val- 
ues of TTiQ, rn-Yji, the squark and gluino masses, and the 
NLO cross section. The quoted systematic uncertainties 
are discussed below. 

The uncertainty coming from the JES corrections is 
one of the most important. It is typically of the order of 
13% for the SM backgrounds and 10% for the signal effi- 
ciencies. The uncertainties on the jet energy resolution, 
on the jet track confirmation, and on the jet reconstruc- 
tion and identification efficiencies were evaluated. They 



lead to systematic uncertainties of 3.5%, 4.0% and 5.4% 
in the "dijet," "3-jets," and "gluino" analyses, respec- 
tively. The trigger was found to be fully efficient for 
the event samples surviving all analysis cuts. Conserva- 
tively, a 2% uncertainty was set on the trigger efficiency. 
The uncertainty on the determination of the luminosity 
is 6.5% [13. All of these uncertainties are fully correlated 
between signal and SM backgrounds. A 15% systematic 
uncertainty was set on the M^/Z-|-jets and tt NLO cross 
sections. The uncertainty on the signal acceptance due 
to the PDF choice was determined to be 6%, using the 
forty-eigenvector basis of the CTEQ6.1M PDF set [l^. 

The signal cross sections are very sensitive to the PDF 
choice and to the renormalization and factorization scale, 
/irf. The nominal NLO cross sections, CTnomi were com- 
puted with the CTEQ6. IM PDF and for /i,f = Q, where 
Q was taken to be equal to rtig for gg production, rnq 
for qq and qq productions, and (m, -I- mg)/2 for qg pro- 
duction. The uncertainty due to the choice of PDF was 
determined using the full set of CTEQ6.1M eigenvec- 
tors, with the individual uncertainties added in quadra- 
ture. The effect on the nominal signal cross sections, 
which varies between 15% and 50%, is dominated by 
the large uncertainty on the gluon distribution at high x. 
The effect of the renormalization and factorization scale 
was studied by calculating the signal cross sections for 
^■ti = Q, fJ-if = Q/2 and fi,-{ = 2 x Q. The factor 



o 
o 



25, Ap=0,tan(P)=3,ti<0 
— NLO cross section 
•— Observed limit 
Expected limit 




240 260 280 300 320 340 360 380 

Squark Mass (GeV) 



:m., Ao=0, tan(P)=3,n<0 
p^ NLO cross section 
•— Observed limit 
A Expected limit 




L \' 


m 


„=500 


, A„=0, lan(p)=3, 


H<0 


\. 







NLO cross section 


r \. 




— •— 


Observed limit 




-^^ 


\ 


--A-- 


Expected iimit 




^^\ 


\ 


k 






: '^^^^S^^N 


V 


k 






'"aTT^. 


-,\ 










■■^ 


\?^^ 




DORunll 


310 pb' 




XJ^ 


^ 



Squark-Gluino Mass (GeV) 



260 280 

Gluino Mass (GeV) 



FIG. 2: For tan (3 = 3, Aq = 0, fi < 0, observed (closed circles) and expected (opened triangles) 95% C.L. upper limits on squark- 
gluino production cross sections combining the analyses for mo = 25 GeV (left), m, — rrig (middle), and mo — 500 GeV (right). 
The nominal production cross sections are also shown, with shaded bands corresponding to the PDF and renormalization-and- 
factorization scale uncertainties. 



two on this scale reduces or increases the nominal sig- 
nal cross sections by 15—20%. The PDF and ^rf effects 
were added in quadrature to compute minimum, (Tmin, 
and maximum, (Tmax, signal cross sections. 

No significant excess of events was observed in the data 
with respect to the SM background expectation in any 
of the three analyses. Therefore, an excluded domain in 
the gluino-squark mass plane was determined as follows. 
The three analyses were run over signal MC samples gen- 
erated in the gluino-squark mass plane to compute signal 
efficiencies. Then, to take advantage of the different fea- 
tures of the three analyses, they were combined in the 
limit computation, with the small overlaps taken into ac- 
count. In the data, no events were selected by more than 
one analysis. 

Limits at the 95% C.L. were computed for three hy- 
potheses on the signal cross sections: nominal, minimum, 
and maximum. Figure |21 shows the observed and ex- 
pected upper limits on squark-gluino production cross 
sections for the three benchmark scenarios. For the "3- 
jets" and "gluino" analyses, the expected limits com- 
puted with the numbers of events reported in Table ^ 
are almost identical to the observed ones. Once the com- 
bination of analyses is performed, the expected limits 
become slightly better than the observed limits at large 
mo and for rriq = nig, as can be seen in Figure [2 

Figure 13 shows the excluded domain in the gluino- 
squark mass plane. The absolute lower limits on the 
squark and gluino masses obtained in the most conserva- 
tive hypothesis, amin, are 325 GeV and 241 GeV, respec- 
tively. The corresponding expected limits are 330 GeV 
and 246 GeV. Table Hill summarizes these absolute limits 
as a function of the signal cross section hypothesis. Lim- 
its were also derived for the particular case rriq = rrig. 
For (Tinin, squark and gluino masses below 337 GeV are 
excluded, while the expected limit is 340 GeV. The ob- 
served limit becomes 351 GeV for CTnom, and 368 GeV for 



These results improve on the previous direct limits on 
squark and gluino masses |3, y, lUl • They were obtained 
within the mSUGRA framework with tan/3 = 3, Aq = 0, 



TABLE III: Absolute lower limits at the 95% C.L. on the 
squark and gluino masses (in GeV) as a function of the choice 
of signal cross section hypothesis as defined in the text. Num- 
bers in parentheses correspond to the expected limits. These 
limits are valid for the mSUGRA parameters: tan/3 — 3, 
Aq = 0, /i < 0. 



Hypothesis 



Gluino mass 



Squark mass 



^nom 



241 (246) 
257 (261) 
274 (280) 



325 (330) 
339 (344) 
352 (358) 



and /i < 0. A general scan of the mSUGRA parameter 
space is beyond the scope of the current analysis, but it 
has been verified that similar results would be obtained 
for a large class of parameter sets. The limits obtained at 
LEP on the chargino (x^) and slepton {£) masses can be 
turned into constraints on the mSUGRA parameters ttiq 
and 1711/2 |l8| - and hence on the squark and gluino masses 
as shown in Fig. |31 The limits from Higg s boson searches 
at LEP are even more constraining |l8| . actually ruling 
out all of the squark and gluino mass domain to which 
the Tevatron could be sensitive. The interpretation of 
these indirect constraints is however more sensitive to 
the details of the model considered than the direct limits 
presented here. 

In summary, a search for events with jets and large 
]^T has been performed in a 310 pb~ data sample from 
pp collisions at 1.96 TeV, collected by the D0 detector. 
Three analyses were designed, specifically targeted to the 
dijet, three-jet, and multijet topologies. The numbers 
of events observed are in agreement with the SM back- 
ground predictions. The results have been interpreted in 
the framework of minimal supergravity with tan/3 = 3, 
Aq = 0, fi < 0. For the central choice of PDF, and 
for a renormalization and factorization scale equal to the 
mass of the squark or gluino produced, the lower limits 
on the squark and gluino masses are 339 and 257 GeV 
at the 95% C.L. Taking into account the PDF uncer- 
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FIG. 3: In the gluino and squark mass plane, excluded regions 
at the 95% C.L. by direct searches in the mSUGRA frame- 
work with tan/3 — 3, Aq = 0, n < 0. The new region excluded 
by this analysis in the most conservative hypothesis (amin) is 
shown in dark shading. The thick line is the limit of the 
excluded region for the CTnom hypothesis. The corresponding 
expected limit is the dotted line. The band delimited by the 
two dashed lines shows the effect of the PDF choice and of a 
variation of /irf by a factor of two. Regions excluded by previ- 
ous experiments are indicated in light shading 0,13, [I3. The 
two thin lines indicate the indirect limits inferred from the 



LEP2 chargino and slepton searches |lq|. The region where 
no mSUGRA solution can be found is shown hatched. 



tainties and allowing for a factor of two in the choice of 
scale, these limits are reduced to 325 and 241 GeV, re- 
spectively. These are the most constraining direct limits 
on the squark and gluino masses to date. 
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